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ABSTRACT 


This  paper  Investigates  the  application  of  the  Kalman  Filter  and  the 
General  Exponential  Smoothing  techniques  of  forecasting.  Both  methods 
are  derived  and  the  similarities  and  differences  between  them  are  dis¬ 
cussed.  The  two  techniques  are  then  applied  to  the  practical  problem  of 
predicting  weekly  losses  suffered  by  the  U.  S.  Marine  Corps  units  in  the 
I  Corps  Tactical  Zone  In  the  Republic  of  Vietnam.  The  mean  absolute 
error  of  the  prediction  is  used  as  the  criterion  for  choosing  the  better 
of  the  two  methods.  Results  are  given  for  both  techniques  as  well  as 
for  the  method  of  linear  regression.  In  general  the  Kalman  Filter  pro¬ 
vides  the  smallest  mean  absolute  error  for  the  three  mathematical  models; 
linear,  growing  sine  with  harmonics  and  frequency  of  sixteen,  thirty-two, 
and  fifty- two  weeks,  and  a  constant  model. 
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I.  INTRODUCTION 


One  of  the  problems  in  any  forecasting  technique  is  what  to  do  about 
additional  data.  That  Is,  given  a  forecast  based  upon  some  data,  how 
does  one  go  about  updating  this  forecast  when  additional  data  is  obtained? 
Two  such  methods  are  discussed  In  some  detail.  These  are  the  methods 
known  as  the  Kalman  Filter,  named  for  R.  E.  Kalman,  and  General  Exponen¬ 
tial  Smoothing,  a  name  coined  by  Robert  G.  Brown. 

Both  the  Kalman  Filter  and  General  Exponential  Smoothing  are  derived 
and  comparisons  made  for  any  similarities  and  differences.  Both  tech¬ 
niques  are  then  applied  to  the  problem  of  forecasting  losses  Incurred  by 
the  U.  S.  Marine  Corps  operating  in  the  I  Corps  Tactical  Zone  of  the 
Republic  of  Vietnam.  The  actual  losses  are  expressed  in  weeks  and  cover 
the  period  from  December  24,  1965,  tc  April  6,  1968.  Losses  Include  all 
Marine  Corps  personnel  who  die  as  a  result  of  battle  as  well  as  non-battle 
causes.  It  also  Includes  those  personnel  who  must  be  evacuated  from 
Vietnam  because  of  wounds  or  injuries  received  from  any  cause. 

The  basis  for  comparing  the  Kalman  Filter  and  General  Exponential 
Smoothing  lies  in  their  similarities.  In  each  case,  a  mathematical  model 
is  assumed.  The  model  consists  of  coefficients,  which  are  unkown,  and 
functions,  called  fitting  functions,  whose  values  are  known.  The  problem 
becomes  one  of  estimating  the  coefficients,  based  upon  previous  data, 
then  using  this  estimate  with  the  fitting  functions  to  obtain  a  forecast. 
The  forecast  error,  the  difference  between  what  actually  occurs  and  what 
has  been  forecasted  to  occur.  Is  weighed  to  obtain  a  new  estimate  of  the 
coefficients. 
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A.  THE  MODELS 

The  models  used  for  forecasting  U.  S.  Marine  losses  are  a  linear 
model , 

C(t)  *  a1  +  agt, 
a  growing  sinusoidal  with  harmonics, 

C(t)  «  (a1  +  a2t)  +  (a3  +  a5t)s1r>^- + 

(a4  +  a6t)cos~-  +  a?  sir4~£  +  a8  cos^-  , 

and  a  constant  model , 

C(t)  ■  a-|  . 

In  the  second  model,  the  forecasts  were  made  for  values  of  w  of  sixteen, 
thirty- two,  and  fifty- two  weeks. 

The  second  model  was  chosen  since  there  appeared  to  be  some  cycle  in 
the  number  of  losses  Incurred.  It  was  felt  that  losses  would  Increase 
as  a  result  of  Increased  enemy  activity.  Such  activity  was  at  a  low  ebb 
during  the  time  of  the  monsoon  rains.  During  the  dry  period,  the  tempo 
of  enemy  offensive  actions  Increased  for  a  time  with  a  resultant  increase 
in  Marine  Corps  losses. 

B.  NOTATION 

The  notation  used  above  reflects  that  used  throughout  this  paper. 

The  symbols  e(t)  or  C(t)  represent  the  forecast  of  losses  for  the  tth 
time  period.  The  same  symbols  without  the  *  above  them  are  the  observed 
value  or  the  actual  losses  during  the  tth  period.  It  is  also  convenient 
to  adopt  vector  and  matrix  notation.  Thus,  an  alphabetical  symbol,  such 
as  g,  will  represent  a  column  vector,  while  aT  will  denote  a  row  vector. 
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In  general,  a  vector  with  a  subscript  will  represent  the  vector  at  a 
specific  point  In  time.  Matrices  will  be  denoted  by  a  capital  alpha¬ 
betical  character.  In  the  event  the  matrix  symbol  Is  subscripted,  it 
will  denote  the  matrix  obtained  by  taking  the  expected  value  of  a  col¬ 
umn  vector,  the  first  subscript,  multiplied  by  the  transpose  of  the  vector 
denoted  by  the  second  subscript.  A  matrix  with  a  single  subscript  Is  the 
expected  value  of  the  vector  represented  by  the  subscript  times  Its  trans¬ 
pose.  For  example, 

Cab  -  E(abT),  and 
Ca  *  E(aJ). 

II.  THE  KALMAN  FILTER 

Basic  to  every  forecasting  technique  is  the  criterion  which  is  used 
to  determine  the  forecast.  In  the  case  of  the  Kalman  Filter,  which  esti¬ 
mates  the  coefficients  of  the  model,  the  criterion  Is  to  minimize  the 

trace  of  the  estimation  error  covariance  matrix,  C„,  where  e  Is  the 

e 

unobservable  error  vector  between  the  true  value  of  the  model  coeffi¬ 
cients  and  the  current  estimate  of  these  coefficients.  This  estimate  is 
the  minimum  mean  square  error  estimator. 

A  fundamental  theorem  of  estimation,  called  the  Gauss-Markov  Theorem, 
is  very  important  In  any  form  of  estimation.  The  theorem  shall  be  stated 
but  no  proof  is  provided. 

Theorem:  If  e  is  the  vector  of  the  observed  data  at  times  t^  and  i 
Is  a  linear  function  of  @,  x  and  g  are  random  variables  with  moment  matrices 
C  ,  C  ,  and  C  ,  and  if  C"^  exists,  then  the  linear  minimum  mean  square 

A  A0 

estimate  x  of  x  given  e  is 


?  ■  v1?  • 

an<*  Ce  “  cx  '  cxeC9  cxe  * 

Throughout  this  paper,  the  data  will  be  assumed  to  be  linearly 
related  to  the  vector  x  (or  a).  Thus, 

e  *  Bx  +  v, 

where  v  is  the  measurement  error  or  noise  vector  and  B  is  an  m  x  n 
matrix  whose  rows  are  the  vectors  fT(t),  the  vector  of  fitting  functions 
evaluated  at  time  t.  For  example,  if  e(t)  *  x1  +  x2t  +  v,  the  vector 
of  fitting  functions  fT(t)  *  (l,t)  and  the  vector  of  coefficients 
xT  -  (x1  ,  x2). 


Then, 


Cxe 

-  E[x(B*+v)T] 

or 

C  Q 

=  c  bt  +  , 

xe 

X  XV 

and 

Ce  = 

E[B^  +  x)(Bx.+  y)T] 

or 

Ce  = 

BC  BT  +  BC  +  CT  B 

X  XV  XV 

From  the  fundamental  theorem, 

x  *  [C  BT  +  C  1[BC  BT  +  BC  +  CT  +  C  I"1*  , 

*  L  x  xvJL  x  xv  xv  vJ  §  * 

and  C  *  C  -  [C  BT  +  C  ][BC  BT  +  CT  BT  +  BC  +  C  I'Ve  BT  +  C  ]T. 

e  x  L  x  xvJL  x  xv  xv  vJ  L  x  xvJ 

Normally  the  vectors  x  and  v  are  assumed  to  be  independent.  This, 
with  the  assumption  that  E(v)  -  0,  Implies  that  Cxy  Is  identically  zero. 
With  this  assumption,  the  above  equations  reduce  to 

x  *  CXBT[BCXBT  +  Cy]e 
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C0  ■  c  -  CvBT[BC  bt  +  cj-Wf  . 


By  the  use  of  the  matrix  Identities 

(2.1)  ABT[C  +  BAB1]"1  -  [A"1  +  BTC"1B]"16TC",1» 
and  (2.2)  [C"1  +  bV1#]’1  -  C  -  CBT[BCBT  ♦  Aj^BC, 


the  following  equations  are  obtained: 

x  -  [cj1  +  B^BrVc^e 

and  Cfi  »  [CJ1  ♦  B^B]*1. 


There  still  remains  the  matrix  Cx  which  must  be  known  before  the 

A 

estimate  x  can  be  made.  In  general  little  Information  Is  available 
concerning  x.  In  this  case  It  Is  reasonable  to  assume  that  the  diagonal 
elements  of  Cx  are  larger  than  any  finite  riL^er.*  With  this  assumption, 
C'1  Is  considered  to  be  0.  Then, 

x  -  Cbtc;1b]“1btc;1©, 

and  Ce  «  [B^B]'1. 

A.  THE  KALMAN  FILTER  EQUATIONS 

The  foregoing  has  been  the  groundwork.  Suppose  that  m  observations, 

A 

9m,  have  been  taken  and  xm  Is  the  minimum  mean  square  error  estimate 
based  upon  the  m  observations.  The  matrix  Is  available.  In  addition, 
suppose  that  r  additional  observations  are  taken.  The  problem  Is  to 


determine  the  estimate  x  . 

A  ^ 

Since  xm  is  available  and  Is 


xm  «  [BmT(c”)’1Bm]"1BnT(c”)"1  ?m 


cJJ  ■  [BmT(cJ)'1Bm]"1  . 

Also,  em  -  Bmxm  +  vm  , 

<v  «• 

and  em+r  ,  Bm+r  xm+r  +  ym+r  . 

The  noise  vector  vm  Is  Independent  of  the  vector  vr.  This  Implies  that 
Cyfr  is  a  diagonal  matrix  with  elements  cj  and  cj.  By  partitioning  the 
matrix  Bm+r  Into  components  Bm  and  Br  and  the  vector  8nHT  Into  em  and  er, 
the  new  estimate 

i"fr  -  [<*''  ♦  B'VrVr’KO'1*" +  BrT(Cv>"1or3 

Is  obtained.  Also, 

cf r  -  [(C^r1  ♦  BrT(cj rVr1 . 

Then  define  K  *  cJfrB,rT(cJr1  .  From  equation  (2.1), 


K  *  C™BrT[cJ  +  BrcjjBrT]”1  . 


From  equation  (2.2) , 


C?r  .  <£-KBr<£. 


Finally,  the  new  estimate  Is 


j^H-r  „  jjm  +  K^0r  _  Br^m-j 


These  last  three  equations  are  the  Kalman  Filter.  [Liebelt  1967, 

pp.  165-166]  Note  that  succeeding  estimates  xnH>r  are  computed  recursively 

from  the  prior  estimate  x1". 
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B.  REMARKS 

Notice  that  with  the  assumption  that  Cy  Is  known,  neither  K  nor 
depend  upon  the  data  e.  These  elements  could  be  precomputed  and  re¬ 
tained  for  later  use  as  needed.  Also,  It  Is  not  necessary  to  assume 
that  the  matrix  C”1  Is  zero.  If  one  has  Information  about  the  distri¬ 
bution  of  x,  the  elements  of  Cx  can  be  computed  and  used  In  determining 

v 

The  vector  K  may  be  thought  of  as  a  weighting  factor  applied  to  the 
forecast  error.  Notice  that  the  weight  for  any  error  Is  proportional 
to  the  Inverse  of  the  variance  of  the  forecast.  A  large  variance  In 
the  measurement  or  noise  for  any  particular  forecast  will  cause  the 
error  to  be  given  less  weight  In  determining  the  estimate  of  the  coef¬ 
ficients  to  be  used  In  computing  the  next  forecast. 

The  dimension  of  the  matrix  Cy  Is  always  r  x  r.  If  measurements 
or  observations  are  taken  at  each  succeeding  time  period,  Cy  Is  a  scalar 
quantity  and  Is  the  variance  of  the  measurement  for  that  period.  In  the 
event  that  the  matrix  Cy  Is  diagonal  with  the  diagonal  elements  equal  to 
some  constant  raised  to  some  power,  the  matrix  Cfi  will  reach  a  steady 
state.  That  Is, 

C"*r  -  C*  . 
e  e 

The  power  to  which  the  constant  Is  raised  corresponds  to  the  time  period 
for  which  the  observation  Is  made. 
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III.  GENERAL  EXPONENTIAL  SMOOTH I N6 


As  Is  the  case  In  the  Kalman  Filter  method*  the  linear  model 

x  ■  Ae 

Is  assumed,  but  rather  than  estimating  x  In  order  to  minimize  the 
diagonal  elements  of  some  matrix.  It  Is  desired  to  minimize  the  sum  of 
the  weighted  squared  error, 

z  w2[e(j)  -  e(d)]2  . 
j-1  3 

If  W  Is  a  diagonal  matrix  with  elements  (w-j , . . .  ,wm) ,  then 

z  w?[e(1)  -  e(1)]2  »  [H(e  -  e)]TW(e  -  0)  . 

1-1  1  '  '  '  - 

Also,  It  Is  assumed  that 

*  a 

0  *  8x  , 

where  the  matrix  B  Is  defined  the  same  as  In  the  discussion  of  the  Kalman 
Filter. 

Define  the  matrix  F  to  be  BTWWTB.  Then 
[W(e  -  0)]TW(e  -  e)  «  [eTAT  -  eTWWTBF-1]F[eTAT  -  oWbf"1]1 

-  0TWTBF_1[eTWwVT  +  eTwTWe  . 

The  matrix  F  Is  positive  definite.  Therefore  the  choice  of  AT  which 
minimizes  the  above  expression  Is 

at  -  hhtbf_1  . 
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Finally, 


xT  -  eWBF-1  * 

Notice  that  this  Is  Identical  with  the  6auss-Martov  estimate  with 
F  -  bV’b  and  WWT  •  C"1 

V  w  • 

A.  THE  MEI6HTIN6  VECTOR 

A  fundamental  hypothesis  of  Exponential  Smoothing  Is  that  the 
vector  of  fitting  functions  f(t+l )  must  be  a  linear  combination  of  f(t). 
The  matrix  of  coefficients  of  this  combination,  which  does  not  depend 
on  time.  Is  called  the  transition  matrix,  L. 


From  this. 

f(t+l) 

-  Lf(t)  . 

and 

e(t+l) 

-  fT(l)LTx(t-l) 

where 

x(t)  * 

(LT)^x(0)  . 

The  fitting  functions  are  always  computed  at  the  number  of  time  periods 
in  the  future  for  which  the  forecast  Is  desired.  This  Implies  that  the 
present  period  Is  always  considered  as  the  period  t  «  0.  Then  the 

M  A 

problem  becomes  how  to  compute  x(t)  recursively  from  x(t-l )  so  that 

«» 

z  b*[e(t-1)  -  xftjff-l)]2  Is  minimized. 

1*0 

Define  the  matrix 

T  1  T 

F(T)  -  z  bVHJf'M)  , 

1-0  ' 

and  the  matrix  F  to  be  the  limit  of  F(t)  as  t  becomes  very  large.  The 
matrix  F  will  exist  If  each  fitting  function  fj(t)  Is  greater  than  the 
square  root  of  b“t. 
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Let  the  data  vector  be 


g(T)  «  i  b  e(T-1)f(-1)  . 

1-0 

Then  the  1th  element  of  the  vector  g(T)  Is  z  bke(T-k)f,(-k) 

k«0  1 


g(T)  -  e(T)f(0)  +  bL’  g(T-l)  . 


F(T)x(T)  -  z  bVc-Df^-Dxd)  . 


e(T-1 )  -  f'(-1)x(T)  . 


Hence, 


F(T)i(T)  -  g(T)  . 


Substituting  the  steady  state  matrix  F  and  combining  the  two  equation 
for  g(T),  the  expression 

x(T)  *  0(T)F'1f(O)  +  bF"1L_1Fx(T-l) 

Is  obtained.  [Brown  1963]. 

But  bF'W  «  LT  -  F_1f(0)fT(0)LT 


fT(1)  «  fT(0)LT  . 


Define  the  vector  h  to  be  F”  f(0).  Then  the  expression  for  x(T)  becomes 


x(T)  «  he(T)  +  [LT  -  hfT(l)]x(T-l)  , 


or  x(T)  «  L'x(T-l)  +  h[e(T)  -  e(T)j  . 

This  last  equation  Is  the  Exponential  Smoothing  means  of  estimating  the 
model  coefficients.  [Brown  1963,  pp.  174>177].  The  vector  h  is  the  weight 
of  the  forecast  error  applied  in  obtaining  the  new  estimate  x(T).  Note 
that  x(T)  Is  confuted  In  a  recursive  way  from  x(T-1). 
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B.  THE  FORECAST  BIAS 

One  My  question  the  eccuracy  of  the  forecast  obtained  by  Mans  of 
Exponential  Snoothlng.  It  has  been  shown  [Bessler  and  Zehna  1968]  that 
for  the  nodel  C(t)  »  a  +  e(t)  , 

where  E[C(t)]  »  a 

that  the  forecast  Is  biased  In  the  early  stages  by  a  factor  of  l-b*, 
where  b  Is  the  weighting  scalar  applied  to  past  data.  Since  b  Is 
between  zero  and  one,  then  E[C(t)]  -  a 

when  bt  has  vanished.  However,  this  shows  that  even  where  the  nethe- 
mtlcal  nodel  Is  correct,  there  Is  an  Initial  bias  In  the  forecast. 

This  nay  lead  one  to  ponder  what  occurs  when  the  assuMd  nodel  Is 
Incorrect.  It  has  been  demonstrated  [Brawn  1963,  p.  128]  that  In  the 
case  where  the  assuned  nodel  Is 

E[C(t)]  -  a 

but  the  data  Is  actually  corresponding  to  the  nodel 
E[C(t)J  »  a1  +  a2t  . 

that  E[C(t)]  ■  a<|  ♦  a2t  -  a^O-h)”1  • 

after  bt  vanishes.  Thus,  even  after  a  sufficient  amount  of  tine  has 
elapsed  so  that  b*  ■  0,  the  forecast  Is  Incorrect  or  biased  by  a  constant 
tines  the  slope  of  the  line  through  the  data. 

The  next  step  was  to  determine  If  these  results  could  be  generalized. 
To  do  this,  the  Exponential  Snoothlng  equation 

i(t)  -  LT»(t-1)  +  h[C(t)  -  fT(l)*{t-l» 

was  examined.  This  can  be  rewritten  as 


a(t)  -  tz1  [LT  -  hfT(l)]khC(t-k)  +  [LT  -  hfT(1  )]*•«)). 
k-0  "  " 

Then  If  E[C(t)3  -  gT(t)a  , 

E[C(t)]  -  fT{l)tz\lT-hfT(1)]llhgT(t-k)a  +  fT(l)[LT-hfT(l)]tE[a(0)]. 
k-0 . 

The  above  expression  Is  true  In  general  when  the  assumed  model  Is 
C(t)  -  fT(t)a  +  e(t) 
and  the  data  actually  follows  the  model 

C(t)  «  gT(t)a  +  e(t). 

.  *  t  -w 

In  order  to  compare  this  expression  with  those  obtained  previously  for 
the  constant  model ,  the  values 

L  -  1, 
f(t)  «  1, 
and  g(t)  -  1 

are  substituted  to  obtain 

.  t~1  L 

E[C(t)3  -  i  [l-h]*h  a  . 

k-0 

This  can  be  written  as 

EtCtt)]  -  , 

where  b  1$  1-h.  Thus  the  result  Is  the  same  as  before. 

Next,  It  was  desired  to  consider  the  situation  where  the  data  was 
modeled  by  E[C(t)]  -  a  ,  but  the  data  actually  obeyed  the  model 
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E[C(t)  *  a-  +  a2t. 


Then  substituting  the  values 
LT  -  1, 
fT0)  -  i. 

9T(t)  *  O.t). 

?T  *  (*V  *2> 


into  the  expression  for  E[C(t)3, 


E[C(t)3  «  i  [1-h]k(l.t)(ai.  aJT  h  . 
k*0  1  z 

This  becomes  E[C(t)3  »  (a-j+agtHl-h*)  -  a2(l-h)lT'(l-h*) 


When  hx  vanishes,  the  expression  becomes 


E[C(t)]  «  a]  +  a2t  -  a2(l-h)h  , 

which  Is  Identical  with  the  prior  result  with  b  •  1-h  . 

Notice  that  the  expression  for  E[C(t)3  has  been  used  to  determine  If 
a  bias  exists  when  the  correct  mathematical  model  Is  assumed  as  well  as 
to  determine  the  mount  of  the  bias  when  an  Incorrect  model  Is  assumed. 

It  was  desirable  to  determine  If  the  forecast  was  biased  when  the  .assumed 
linear  model  was  correct.  The  equation  then  becomes 

E[C(t)3  -  (a1+a2t)fT(l)[I-LT+hfT(l)3"1[I-aT-hfT(l))t3h 

-  «2fT(l  )D-lT+hfT(l )Y}  [lT-hfT{l  )3Cl-LT+hfT(l  J3’1 

Cl-(LT-hfT(l))t"13h-(t-l)n-LT+hfT(l)3*,[LT-hfT(l)3t*1h  . 
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■  -  III... I  . .  . . .  M  - - - — 


It  became  Important  to  know  If  the  vector  [LT-hfT(l)]th  vanished  as 
t  became  large.  If  the  vector  falls  to  vanish,  the  forecast  will  remain 
biased.  The  vector  will  vanish  If  the  eigenvalues  of  [LT»hfT(l)]  are 
all  between  minus  one  and  plus  one.  To  demonstrate  this,  let 
A  *  LT-hfT(l)»  where  A  Is  n  x  n  and  nonsingular.  Also,  let  the  eigen¬ 
values  of  A  be  m^,  mg •  •  •  •*R|n*  where  each  m..  lies  between  minus  one 
and  plus  one.  Let  v^,  v2,  .  .  .,vn  be  eigenvectors  of  A.  These  vectors 
form  a  basis  for  En. 

Then  Av1  »  m^  , 
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In  the  case  of  the  linear  model,  the  eigenvalues  are  equal  to 

* 

1  -  (hj  +  h2)2’1.  These  eigenvalues  Mill  be  between  minus  one  and  plus 
one  If  hj  +  h2  lies  between  zero  and  four.  The  h  for  the  linear  model 
Is  such  that  this  Is  true,  If  b  Is  between  zero  and  one.  To  demonstrate 
this,  recall  that 

h  -  F“1f <0) . 

F  -  e  bVf-Of'H)  , 

1-0  - 


and 

fT(t)  - 

0.  t)  . 

Then 

F-'  - 

Tv 

(1-b)2 

a 

(1-b)2 

(l-b)V1 

Then  hT  *  (h-j ,  h2) 


[1-b2,  (1-b)2].  Hence, 


'  h,  +  h2  »  2(l-b)2, 

which  Is  les's  than  four  and  Is  greater  than  zero. 

Therefore,  the  vector  [lT-hfT(l )]*h  vanishes  In  the  linear  model  for 
all  values  of  b  which  are  between  zero  and  one.  Then 

E[C(t)]  =  (a^+a2t)fT(l  )[I-LT+hfT(l  )]*,1h 

-  agf^DCl-J+hf^Dr^J-hf^DlCl-J+hf^Dr1  h  . 

Thus  the  forecast  Is  unbiased  If 

fT(l)[I  -  LT  +  hfT(l )3_1  h  -  1 

and  fT(l)[I-LT+hfT(l)]‘1[LT-hfT(l)][I-LT+hfT(l)r1  h  -  0  . 

For  a  specific  example  see  Chapter  V. 
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C.  EXPONENTIAL  SMOOTHING  COMPARED  WITH  KALMAN  FILTER 

The  Exponential  Smoothing  equation  is  a  special  case  of  the  Kalman 
Filter  In  the  steady  state.  Recall  that  the  matrix  F  was  defined  to  be 
BTWWTB.  If  WWT  is  equivalent  to  the  matrix  C”1  used  in  the  Kalman 
Filter  equations,  then  F"^  is  Identical  to  the  matrix  Cg  derived  in 
Chapter  II.  Thus,  In  the  case  of  fitting  functions  consisting  of  poly¬ 
nomials,  exponentials,  and  sinusoldals,  which  can  be  generated  by  a 
transition  matrix,  the  Kalman  Filter  estimate  becomes  identical  with 
that  ot  Exponential  Smoothing  when  C  reaches  the  steady  state.  However, 

C 

since  the  Kalman  Filter  method  does  not  require  that  the  fitting  func¬ 
tions  be  generated  by  means  of  a  transition  matrix,  it  may  be  applied  to 
models  other  than  those  listed  above. 

The  primary  difference  between  the  Kalman  Filter  and  Exponential 
Smoothing  Is  the  models  for  which  they  are  valid.  The  Kalman  Filter  is 
valid  both  In  a  dynamic  system  as  well  as  in  the  steady  state.  The 
Exponential  Smoothing  technique  Is  strictly  valid  only  in  the  steady 
state. 
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IV.  FORECASTING  MARINE  LOSSES 

In  order  to  see  how  the  Kalman  Filter  and  the  Exponential  Smoothing 
techniques  compared,  both  methods,  along  wltn  linear  regression,  were 
used  to  forecast  0.  S.  Marine  Corps  losses  in  the  Republic  of  Vietnam. 

The  actual  data  used  Included  one  hundred  and  nineteen  weeks  from 
December  26,  1965,  through  April  6,  1968.  The  losses  Included  those 
Marines  who  either  were  killed  or  were  wounded  and  evacuated  from 
Vietnam. 

Three  models  were  Investigated.  These  were  a  growing  sinusoidal 
with  harmonics  with  a  cycle  of  sixteen,  thirty-two,  and  fifty-two  weeks, 

C(t)  =  a.j+a2t+(a3+a5t)5ln(pt)  +  (a4+agt)cos{pt) 

+  a7sin(2pt)  +  a8cos(2pt)  , 

a  linear  model , 

C(t)  =  a^  +  a2t  , 

and  a  constant  model , 

C(t)  *  a  . 

The  objective  was  to  minimize  the  mean  absolute  error  where  the  absolute 

a  l  a 

error  for  the  t  forecast  was  the  absolute  value  of  C(t)  -  C(t).  The 
mean  error  and  the  estimated  variance  of  the  errors  were  Included  for 
comparison. 

Both  the  Kalman  Filter  and  Exponential  Smoothing  require  an  estimate 
of  the  coefficients  a(0)  in  order  to  begin  forecasting.  Captain  Paul 
William  O'Brien,  U.  S.  Marine  Corps,  had  studied  the  problem  of 
Exponential  Smoothing  In  forecasting  Marine  monthly  losses.  [O'Brien  1968]. 


23 


Using  a  least  squares  line  through  the  monthly  losses  from  March  1965 
to  December  1966  and  the  least  squares  Fourier  fit  through  the  deviations 
about  the  regression  line.  Captain  O'Brien  obtained  Initial  values  of  the 
coefficients  for  the  growing  sinusoidal  model.  The  coefficients  of  the 
least  squares  line  through  the  data  were  used  as  the  initial  values  in  the 
linear  model.  These  values  were  scaled  by  one-fourth  and  used  as  the  Ini¬ 
tial  values  for  forecasting  weekly  losses.  Thus  the  growing  sinusoidal 
model  was  Initialized  with 

aT(0)  =  (-36.45,  15.675,  62.61,  31.3975,  .6325,  .6325, 

-35.1175,  54.76)  , 

and  the  linear  model  with 

aT(0)  -  (-36.45,  15.675)  . 

A.  FORECASTING  WITH  EXPONENTIAL  SMOOTHING 

In  each  model  the  value  of  the  weighting  factor  by  which  the  errors 
were  discounted  was  chosen  to  be  the  solution  to 

1  -  bn  =  .25  , 

where  n  was  the  number  of  coefficients  in  the  model.  The  weighting 
vector  h  used  In  the  Exponential  Smoothing  equation  was  computed  by 
determining  the  steady  state  matrix  F  [Brown  1963].  Then 

h  *  F'1  f(0)  . 

The  transition  matrix  In  each  model  was  determined  by  Inspection, 
remembering  trlgnometrlc  Identities  for  s1n(A+B)  and  cos(A+B). 

The  forecasting  procedure  Involved  several  steps.  These  were: 

Step  li  Initialize  the  vector  a(0)  and  compute  f(1). 


24 


Step  2;  Compute  the  forecast  C(l)  by 
C(1)  »  fTO)a(0). 

Step  3;  Compute  the  forecast  error 
e(l)  «  C(!)  -  CO). 

Step  4:  Obtain  a  new  estimate  a(l)  by 
a(1)  *  LTa(0)  ♦  he(l). 

Step  5;  Compute  the  forecast  C(2)  by 
C(2)  -  fT(l)aO). 

Step  6;  Repeat  Steps  3  through  5  for  each  successive  forecast 
forecast  with  the  appropriate  time  arguments. 

The  first  model  investigated  was  the  growing  sinusoidal  with 


harmonics. 

For  this  model  the 

transition  matrix  is 

1 

0 

0 

0 

0  0  0 

0 

1 

1 

0 

0 

0  0  0 

0 

0 

0 

cos  p 

sin  p 

0  0  0 

0 

L  * 

0 

0 

-sin  p 

cos  p 

0  0  0 

0 

0 

0 

cos  p 

sin  p 

cos  p  sin  p  0 

0 

0 

0 

-sin  p 

cos  p 

-sin  p  cos  p  0 

0 

0 

0 

0 

0 

0  0  cos  2p  sin  2p 

0 

_ 

0 

0 

0 

0  0  -sin  2p  cos  2p 

where  p  is  ?w  divided  by  the  cycle  length.  Also 

fT(l)  *  (1,  1,  sin  p,  cos  p,  sin  p,  cos  p.  sin  2p,  cos  2p). 

The  values  of  h  used  were,  for  the  sixteen  week  cycle 
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hT  =  (0.064709,  0.0011416,  0.025939,  0.12583,  0.00046034, 
0.0022276,  0.020951,  0.059457), 

for  a  thirty-two  week  cycle 

hT  *  (0.072061,0.0012067,0.056675,0.12816,0.00096984, 
0.0021717,  0.042191,  0.049775), 

and  for  a  fifty-two  week  cycle 

hT  *  (0.089418,  0.0013569,  0.10333,  0.13127,  0.0016449, 
0.0020314,  0.063796,  0.029308). 

The  results  obtained  for  this  model  for  the  three  cycles  are  contained  in 
Tables  I  through  III. 

The  second  model  was  the  linear  model.  The  same  procedure  was  followed 
as  outlined  above.  The  transition  matrix  was 


and  fT(l)  *  (1,1), 

hT  =  (0.25,  0.01795)  . 

The  results  from  this  model  are  contained  in  Table  IV. 

The  last  model  to  which  Exponential  Smoothing  was  applied  was  the 
constant  model .  With  this  model 

L  -  1, 

fT(1)  *  1, 

and  h  ■  0.25  . 
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TABLE  I 


GENERAL  EXPONENTIAL  SMOOTHING  RESULTS 
GROWING  SINE  WITH  HARMONICS 
CYCLE  LENGTH  16  WEEKS 


WEEK 


LOSSES 


ESTIMATED  LOSSES 


ERROR 


1 

61 

47 

14 

2 

64 

32 

32 

3 

84 

30 

54 

4 

53 

58 

-5 

5 

87 

89 

-2 

6 

116 

120 

-4 

7 

119 

128 

-9 

8 

108 

102 

6 

9 

167 

55 

112 

10 

405 

36 

369 

11 

89 

107 

-18 

12 

85 

107 

-22 

13 

339 

139 

200 

14 

208 

251 

-43 

15 

80 

294 

-214 

16 

116 

270 

-154 

17 

141 

245 

-104 

18 

133 

229 

-96 

19 

125 

227 

-102 

20 

158 

244 

-86 

21 

164 

282 

-118 

22 

82 

319 

-237 

23 

203 

318 

-115 

24 

139 

337 

-198 

25 

133 

321 

-188 

26 

198 

304 

-106 

27 

151 

315 

-164 

28 

140 

322 

-182 

29 

175 

332 

-157 

30 

359 

352 

7 

31 

280 

405 

-125 

32 

1253 

404 

849 

33 

273 

627 

-354 

34 

150 

491 

-341 

35 

257 

348 

-91 

36 

559 

307 

252 

37 

258 

402 

-144 

38 

290 

411 

-121 

39 

382 

418 

-36 

40 

547 

426 

121 

41 

119 

451 

-332 

42 

125 

353 

-228 

43 

85 

317 

-232 

44 

54 

341 

-287 

45 

343 

408 

-65 
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TABLE  I  (CONTINUED) 


WEEK 

LOSSES 

ESTIMATED  LOSSES 

ERROR 

46 

118 

556 

-438 

47 

% 

580 

-484 

48 

80 

541 

-461 

49 

262 

467 

-205 

50 

131 

447 

-316 

51 

85 

408 

-323 

52 

177 

395 

-218 

53 

106 

434 

-328 

54 

435 

444 

-9 

55 

145 

512 

-367 

56 

91 

439 

-348 

57 

126 

338 

-212 

58 

442 

274 

168 

59 

86 

337 

-251 

60 

121 

320 

-199 

61 

153 

340 

-187 

62 

327 

373 

-46 

63 

191 

429 

-238 

64 

220 

403 

-183 

65 

336 

367 

-31 

66 

280 

364 

-84 

67 

237 

357 

-120 

68 

202 

359 

-157 

69 

300 

369 

-69 

70 

273 

407 

-134 

71 

389 

412 

-23 

72 

646 

419 

227 

73 

732 

460 

272 

74 

359 

486 

-127 

75 

400 

397 

3 

76 

253 

362 

-109 

77 

247 

336 

-89 

78 

249 

350 

-101 

79 

248 

373 

-125 

80 

536 

376 

160 

81 

249 

426 

-177 

82 

256 

363 

-107 

83 

254 

324 

-70 

84 

315 

333 

-18 

85 

208 

406 

-198 

86 

269 

464 

-195 

87 

203 

525 

-322 

88 

290 

521 

-231 

89 

608 

495 

113 

90 

357 

517 

-160 

91 

260 

439 

-179 

92 

189 

360 

-171 

93 

166 

317 

-151 

94 

262 

321 

-59 

95 

216 

371 

-155 
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TABLE  I  (CONTINUED) 


WEEK 

LOSSES 

96 

224 

97 

304 

98 

493 

99 

234 

100 

262 

101 

466 

102 

367 

103 

357 

104 

414 

105 

393 

106 

243 

107 

438 

108 

208 

109 

404 

110 

948 

111 

735 

112 

1056 

113 

787 

114 

1339 

115 

453 

116 

563 

117 

475 

118 

676 

119 

582 

Mean  Absolute  Error  * 

170 

Mean  Error  = 

-68.97 

Error  Variance  * 

4.42  X 

ESTIMATED  LOSSES  ERROR 


385 

-161 

366 

-62 

343 

150 

364 

-130 

323 

-61 

338 

128 

440 

-73 

500 

-143 

520 

-106 

505 

-112 

444 

-201 

337 

101 

319 

-111 

277 

127 

333 

615 

529 

206 

585 

471 

656 

131 

590 

749 

667 

-214 

499 

64 

449 

26 

442 

234 

517 

65 
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TABLE  II 


GENERAL  EXPONENTIAL  SMOOTHING  RESULTS 
GROWING  SINE  WITH  HARMONICS 
CYCLE  LENGTH  32  WEEKS 


WEEK 

LOSSES 

ESTIMATED  LOSSES 

ERROR 

1 

61 

60 

1 

2 

64 

64 

0 

3 

84 

63 

21 

4 

53 

67 

-14 

5 

87 

65 

22 

6 

116 

77 

39 

7 

119 

101 

18 

8 

108 

125 

-17 

9 

167 

145 

22 

10 

405 

177 

228 

11 

89 

266 

-177 

12 

85 

241 

-156 

13 

339 

210 

129 

14 

208 

246 

-38 

15 

80 

231 

-151 

16 

116 

176 

-60 

17 

141 

140 

1 

18 

133 

120 

13 

19 

125 

108 

17 

20 

158 

106 

52 

21 

164 

125 

39 

22 

82 

153 

-71 

23 

203 

162 

41 

24 

139 

212 

-73 

25 

133 

240 

-107 

26 

198 

264 

-66 

27 

151 

301 

-150 

28 

140 

316 

-176 

29 

175 

321 

-146 

30 

359 

33^ 

28 

31 

280 

386 

-106 

32 

1253 

404 

849 

33 

273 

681 

-408 

34 

150 

619 

-469 

35 

257 

521 

-264 

36 

559 

466 

93 

37 

258 

505 

-247 

38 

290 

457 

■'167 

39 

382 

433 

-51 

40 

547 

446 

101 

41 

119 

509 

-390 

42 

125 

444 

-319 

43 

85 

395 

-310 

44 

54 

346 

-292 

45 

343 

302 

41 

30 


SiSfSSgggSSggggjgwg 


WEEK 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 


95 


TABLE  II  (CONTINUED) 

ESTimTEDLOSSES 


LOSSES 

118 

96 

80 

262 

131 

85 
177 
106 
435 
145 

91 

126 

442 

86 
121 
153 
327 
191 
220 
336 
280 
237 
20  2 
300 
273 
389 
646 
732 
359 
400 

253 

247 
249 

248 
536 

249 
256 

254 
315 
208 
269 
203 
290 
608 
357 
260 
189 
166 
262 
216 


348 

324 

297 

270 
300 
290 

271 
287 
287 
385 
389 
381 
388 

485 
464 
454 
451 
492 
480 
471 

486 
473 
442 
401 

390 

373 

391 
476 
564 
524 
496 
429 

374 
334 
310 
378 
358 
347 
340 

352 
331 
331 
312 

320 
415 
420 
394 

353 
317 

321 


ERROR 


-230 

-228 

-217 

-8 

-169 

-205 

-94 

-181 

148 

-240 

-298 

-255 

54 

-399 

-343 

-301 

-124 

-301 

-260 

-135 

-206 

-236 

-240 

-101 

-117 

16 

255 

256 
-205 
-124 
-243 
-182 
-125 

-86 

226 

-129 

-102 

-93 

-25 

-144 

-62 

-128 

-22 

288 

-58 

-160 

-205 

-187 

-55 

-105 


31 


TABLE  II  (CONTINUED) 


WEEK 

LOSSES 

ESTIMATED  LOSSES 

ERROR 

96 

224 

320 

-96 

97 

304 

329 

-25 

98 

493 

366 

127 

99 

234 

451 

-217 

100 

262 

446 

-184 

101 

466 

445 

21 

102 

367 

495 

-128 

103 

357 

498 

-141 

104 

414 

489 

-75 

105 

393 

488 

-95 

106 

243 

474 

-231 

107 

438 

417 

21 

108 

208 

423 

-215 

109 

404 

366 

38 

110 

948 

379 

569 

111 

735 

544 

191 

112 

1056 

617 

439 

113 

787 

754 

33 

114 

1339 

776 

563 

115 

453 

926 

-473 

116 

563 

781 

-218 

117 

475 

678 

-203 

118 

676 

564 

112 

119 

582 

531 

51 

Mean  Absolute  Error  *  161 


Mean  Error 
Error  Variance 


-74.73 
3.92  X  104 
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TABLE  III 


WEEK 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 


GENERAL  EXPONENTIAL  SMOOTHING  RESULTS 
GROWING  SINE  WITH  HARMONICS 
CYCLE  LENGTH  52  WEEKS 


LOSSES  ESTIMATED  LOSSES  error 


61 

64 

84 
53 
87 

116 

119 

108 

167 

405 

89 

85 
339 
208 

80 

116 

141 

133 

125 

158 

164 

82 

203 

139 
133 
198 
151 

140 
175 
359 
280 

1253 

273 

150 

257 
559 

258 
290 
382 
547 
119 
125 

85 

54 

343 


63 

73 

79 

87 
83 

88 
102 
114 
123 
147 
235 

219 
203 
260 
269 
238 

220 
211 
200 
187 
185 
185 
161 

175 
169 
163 
177 

176 
174 
183 

244 
273 
569 
535 
466 
435 
491 
451 
426 
432 
484 
408 
349 
294 

245 


-2 

-9 

5 

-34 
4 
28 
17 
-6 
44 
258 
-146 
-134 
136 
-52 
-189 
-1 22 
-79 
-78 
-75 
-29 
-21 
-103 
42 
-36 
-36 

35 
-26 
-36 

1 

176 

36 
980 

-296 

-385 

-209 

124 

-233 

-161 

-44 

115 

-365 

-283 

-264 

-240 

98 
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TABLE  III  (CONTINUED) 


WEEK 

LOSSES 

ESTIMATED  LOSSES 

ERROR 

46 

118 

292 

-174 

47 

96 

275 

-179 

48 

80 

259 

-179 

49 

262 

247 

15 

50 

131 

292 

-161 

51 

85 

296 

-211 

52 

177 

287 

-110 

53 

106 

307 

-201 

54 

435 

302 

133 

55 

145 

389 

-244 

56 

91 

376 

-285 

57 

126 

346 

-220 

58 

442 

327 

115 

59 

86 

397 

-311 

60 

121 

352 

-231 

61 

153 

320 

-167 

62 

327 

302 

25 

63 

191 

335 

-144 

64 

220 

324 

-104 

65 

336 

324 

12 

66 

280 

358 

-78 

67 

237 

372 

-135 

68 

202 

373 

-171 

69 

300 

366 

-66 

70 

273 

389 

-116 

71 

389 

403 

-14 

72 

646 

449 

197 

73 

732 

559 

173 

74 

359 

672 

-313 

75 

400 

657 

-257 

76 

253 

649 

-396 

77 

247 

596 

-349 

78 

249 

545 

-296 

79 

248 

501 

-253 

80 

536 

462 

74 

81 

249 

509 

-260 

82 

256 

466 

-210 

83 

254 

431 

-177 

84 

315 

402 

-87 

85 

208 

394 

-186 

86 

269 

359 

-90 

87 

203 

347 

-144 

88 

290 

320 

-30 

89 

608 

323 

285 

90 

357 

416 

-59 

91 

260 

422 

-162 

92 

189 

398 

-209 

93 

166 

356 

-190 

94 

262 

313 

-51 

95 

216 

306 

-90 
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TABLE  III  (CONTINUED) 


WEEK 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 
119 


LOSSES  ESTIMATED  LOS SE S 


224 

287 

304 

275 

493 

289 

234 

355 

262 

338 

466 

331 

367 

383 

357 

398 

414 

407 

393 

429 

243 

440 

438 

406 

208 

432 

404 

389 

948 

407 

735 

576 

1056 

653 

787 

801 

1339 

840 

453 

1016 

563 

901 

475 

825 

676 

730 

582 

702 

Mean  Absolute  Error  *  154 
Mean  Error  =  -67.78 

Error  Variance  *  3.94  X  10* 


ERROR 

-63 

29 

204 

-121 

-76 

135 

-16 

-41 

7 

-36 

-197 

32 

-224 

15 

541 

159 

403 

-14 

499 

-563 

-338 

-350 

-54 

-120 


TABLE  IV 


GENERAL  EXPONENTIAL  SMOOTHING  RESULTS 
LINEAR  MODEL 


WEEK 

LOSSES 

ESTIMATED  LOSSES 

ERROR 

1 

61 

11 

50 

2 

64 

42 

22 

3 

84 

65 

19 

4 

53 

88 

-35 

5 

87 

94 

-7 

6 

116 

109 

7 

7 

119 

128 

-9 

8 

108 

142 

-34 

9 

167 

148 

19 

TO 

405 

170 

235 

11 

89 

257 

-168 

12 

85 

227 

-142 

13 

339 

201 

138 

14 

208 

258 

-50 

15 

80 

260 

-180 

16 

116 

222 

-106 

17 

141 

203 

-62 

18 

133 

196 

-63 

19 

125 

187 

-62 

20 

158 

177 

-19 

21 

164 

179 

-15 

22 

82 

182 

-100 

23 

203 

159 

44 

24 

139 

178 

-39 

25 

133 

173 

-40 

26 

198 

166 

32 

27 

151 

181 

-30 

28 

140 

177 

-37 

29 

175 

171 

4 

30 

359 

177 

182 

31 

280 

236 

4 

32 

1253 

257 

996 

33 

273 

568 

-295 

34 

150 

504 

-355 

35 

257 

418 

-161 

36 

559 

384 

175 

37 

258 

449 

-191 

38 

290 

406 

-116 

39 

382 

382 

0 

40 

547 

391 

156 

41 

119 

448 

-329 

42 

125 

360 

-235 

43 

85 

295 

-210 

44 

54 

233 

-179 

45 

343 

177 

166 

46 

118 

223 

-105 

36 


WEEK 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 
S3 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 


TABLE  IV  (CONTINUED) 


LOSSES  ESTIMATED  LOSSES 


96 

8C 

262 

131 

85 
177 
106 
435 
145 

91 

126 

442 

86 
121 
153 
327 
191 
220 
336 
280 
237 
202 
300 
273 
389 
646 
732 
359 
400 

253 

247 
249 

248 
536 

249 
256 

254 
315 
208 
269 
203 
290 
608 
357 
260 
189 
166 
262 
216 
224 


189 
157 
128 
162 
148 
124 
134 
120 
210 

190 
159 
146 
232 
190 
168 
161 
209 
204 
209 
248 
261 
257 
243 
262 
269 
308 
416 
524 
491 
477 
422 
377 
343 
317 
385 
349 
324 
304 
307 
277 
273 
250 
259 
363 
365 
337 
294 
254 
254 
239 


ERROR 


-93 

-77 

134 

-31 

-63 

53 

-28 

315 

-65 

-99 

-33 

297 

-146 

-69 

-15 

166 

-18 

16 

127 

32 

-24 

-55 

57 

11 

120 

338 

316 

-165 

-91 

-224 

-175 

-128 

-95 

219 

-136 

-93 

-70 

11 

-99 

-8 

-70 

40 

349 

-6 

-105 

-148 

-128 

8 

-38 

-15 
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TABLE  IV  (CONTINUED) 


WEEK 

LOSSES 

ESTIMATED  LOSSES 

ERROR 

97 

304 

231 

73 

98 

493 

249 

244 

99 

234 

321 

-87 

100 

262 

296 

-34 

101 

466 

286 

180 

102 

367 

340 

27 

103 

357 

351 

6 

104 

414 

357 

57 

105 

393 

377 

16 

106 

243 

387 

-144 

107 

438 

348 

90 

108 

208 

377 

-169 

109 

404 

330 

74 

110 

948 

353 

595 

111 

735 

536 

199 

112 

1056 

609 

447 

113 

787 

761 

26 

114 

1339 

794 

545 

115 

453 

984 

-531 

116 

563 

857 

-294 

117 

475 

793 

-318 

118 

676 

716 

-40 

119 

582 

718 

-136 

Mean  Absolute  Error  =  128 
Mean  Error  x  -2.21 

4 

Error  Variance  *  3.67  X  10 
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TABLE  V 


GENERAL  EXPONENTIAL  SMOOTHING  RESULTS 
CONSTANT  MODEL 


WEEK 


LOSSES  ESTIMATED  LOSSES  ERROR 


1 

61 

2 

64 

3 

84 

4 

53 

5 

87 

6 

116 

7 

119 

8 

108 

9 

167 

10 

405 

11 

89 

12 

85 

13 

339 

14 

208 

15 

80 

16 

116 

17 

141 

18 

133 

19 

125 

20 

158 

21 

164 

22 

82 

23 

203 

24 

139 

25 

133 

26 

198 

27 

151 

28 

140 

29 

175 

30 

359 

31 

280 

32 

1253 

33 

273 

34 

150 

35 

257 

36 

559 

37 

258 

38 

290 

39 

382 

40 

547 

41 

119 

42 

125 

43 

85 

44 

54 

45 

343 

46 

118 

0 

35 

47 

61 

58 

67 

81 

92 

96 

115 

191 

164 

144 
194 
197 
168 

155 
151 

147 
141 

145 

150 
133 

151 

148 
144 
158 

156 

152 
158 
208 
226 
483 
430 
360 
334 
391 
357 
341 
351 
400 
330 
279 
230 
186 
225 


61 

29 

37 
—8 
29 
49 

38 
16 
71 

290 

-102 

-79 

195 

14 

-117 

-52 

-14 

-18 

-22 

17 

19 

-68 

70 

-12 

-15 

54 

-7 

-16 

23 

201 

72 

1027 

-210 

-280 

-103 

225 

-133 

-67 

41 

196 

-281 

-205 

-194 

-176 

157 

-107 
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TABLE  V  (CONTINUED) 


MEEK 

LOSSES 

ESTIMATED  LOSSES 

ERROR 

47 

96 

199 

-103 

48 

80 

173 

-93 

49 

262 

150 

112 

50 

131 

178 

-47 

51 

85 

166 

-81 

52 

177 

146 

31 

53 

106 

154 

-48 

54 

435 

142 

293 

55 

145 

215 

-70 

56 

91 

198 

-107 

57 

126 

171 

-45 

58 

442 

160 

282 

59 

86 

230 

-144 

60 

121 

194 

-73 

61 

153 

176 

-23 

62 

327 

170 

157 

63 

191 

209 

-18 

64 

220 

205 

15 

65 

336 

209 

127 

66 

280 

240 

40 

67 

237 

250 

-13 

68 

202 

247 

-45 

69 

300 

236 

64 

70 

273 

252 

21 

71 

389 

257 

132 

72 

646 

290 

356 

73 

732 

379 

353 

74 

359 

467 

-108 

75 

400 

440 

-40 

76 

253 

430 

-177 

77 

247 

386 

-139 

78 

249 

351 

-102 

79 

248 

326 

-78 

80 

536 

306 

230 

81 

249 

364 

-115 

82 

256 

335 

-79 

83 

254 

315 

-61 

84 

315 

300 

15 

85 

208 

304 

-96 

86 

269 

280 

-11 

87 

203 

277 

-74 

88 

290 

259 

31 

89 

608 

266 

342 

90 

357 

352 

5 

91 

260 

353 

-93 

92 

189 

330 

-141 

93 

166 

295 

-129 

94 

262 

262 

0 

95 

216 

262 

-46 

96 

224 

251 

-27 
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MEEK 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 
119 


TABLE  V  (CONTINUED) 

losses  ESTIMATED  LOSSES  error 


304 

244 

493 

259 

234 

318 

262 

297 

466 

288 

367 

332 

357 

341 

414 

345 

393 

362 

243 

370 

438 

338 

208 

363 

404 

324 

948 

344 

735 

495 

1056 

555 

787 

680 

1339 

707 

453 

865 

563 

762 

475 

712 

676 

653 

582 

659 

Mean  Absolute  Error  *  122 
Mean  Error  »  20.26 

Error  Variance  *  3.48  X  10* 
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In  this  nodel  the  Initial  estimate  of  a(O)  was  taken  to  he  zero.  Also, 

f 

since  the  Initial  forecasts  are  biased  by  1-h  ,  the  method  of  finite 
exponential  smoothing  [Bessler  and  Zehna  1968]  was  used.  With  this 
procedure,  the  weight  given  the  forecast  error  In  determining  the  new 
value  of  a(t)  Is  h(l-h*)”*  rather  than  the  constant  h.  This  method 
Increased  the  initial  forecasts  and  as  h  vanished,  the  forecasts 
became  identical  with  those  using  the  constant  h  throughout.  These 
results  are  In  Table  V. 

B.  FORECASTING  WITH  THE  KALMAN  FILTER 

The  Kalman  Filter  requires,  in  addition  to  the  Initial  value  of 
a(0),  the  value  of  C"1  in  order  to  forecast.  In  the  models  investigated, 
the  matrix  Cy  was  not  known.  However,  since  the  basic  model  was 

C  *  Ba  +  v 

and  in  forecasting  one  time  period  in  advance,  Cy  is  a  scalar,  it  was 
decided  that  Cy  could  be  estimated  as  the  variance  of  the  forecast  error. 
Accordingly,  the  current  estimate  of  the  error  variance  was  used  as  the 
value  of  Cy  in  the  Kalman  Filter  equations. 

The  Initial  value  of  Ce  was  obtained  from  the  equation 

ce  «  [BTCy1B]”'1 . 

This  required  at  least  two  forecasts  In  order  to  estimate  Cy.  So  that 
a  better  estimate  could  be  made,  Cg  was  Initialized  after  ten  time 
periods.  The  forecasts  for  these  ten  periods  were  obtained  by  using 
the  Initialized  value  of  a(0). 

The  forecasting  procedure  was: 

Step  1;  Initialize  a(10)  *  a(0). 
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Step  2;  Obtain  the  forecasts  for  the  first  ten  periods  from 
C(t)  «  bT(t)a(10). 

Step  3;  Compute  the  forecast  errors,  e(t),  for  these  ten 
periods  and  the  mean  error,  e. 

Step  4;  Obtain  an  estimate  of  Cy  from 

i  10  , 

cv  *  TTr  z  M*)  “  • 

v  10  k«l 

Step  5;  Initialize  Ce  by 

ce(io)  »  [btc~^b]-1 . 

Step  6;  Obtain  new  forecast  C(ll)  by 
C(ll)  =  bT(11)a(10). 

Step  7;  Compute  forecast  error  e(ll)  and  the  new  error  mean,  e. 
Step  8;  Obtain  new  estimate  of  Cy  from 

i  11  ? 

Cy  *  e  [e(k)  -  ir  . 

k«l 

Step  9;  Compute  weighting  vector  K  by 

K  -  Ce(10)b(11 )[Cy  +  bT(11)Ce(10)b(ll )]”^. 

Step  10;  Obtain  new  estimate  a(ll)  from 
a(ll)  ■  a( 10)  ♦  Ke(ll). 

Step  11;  Compute  new  value  of  Ce(ll)  from 

Ce(ll)  -  Ce(10)  -  KbT(ll)Ce(10). 
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Step  12;  Repeat  Steps  6  through  11  with  the  appropriate 

time  arguments  for  each  successive  forecast. 

In  general  the  value  of  Cy  will  be  different  from  Cft.  An  exception 

Is  In  the  constant  model  where  they  are  identical.  In  this  model  the 

computation  of  K  In  Step  9  above  involved  the  old  estimate  of  the  error 

variance  as  the  value  of  CA  and  C„  was  the  most  recent  estimate  of  the 

e  v 

forecast  error  variance. 

C.  FORECASTING  WITH  LINEAR  REGRESSION 

The  method  of  Least  Squares  was  applied  to  the  data  using  the 
linear  model.  The  procedure  employed  was  to  fit  the  least  squares  line 
through  the  first  twelve  weekly  loss  figures.  The  coefficients  of  this 
line  were  used  to  obtain  a  forecast  of  losses  for  the  thirteenth  period. 
A  new  line  was  fitted  to  the  thirteen  observations  and  the  new  coeffi¬ 
cients  used  for  the  next  forecast.  This  procedure  was  continued  through 
each  succeeding  time  period. 
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TABLE  VI 


KEEK 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 


KALMAN  FILTER  RESULTS 
GROWING  SINE  WITH  HARMONICS 
CYCLE  LENGTH  16  WEEKS 


LOSSES  ESTIMATED  LOSSES 


61 

64 

84 
53 
87 

116 

119 

108 

167 

405 

89 

85 
339 
208 

80 

ns 

141 

133 

125 

158 

164 

82 

203 

139 
133 
198 
151 

140 
175 
359 
280 

1253 

273 

150 

257 
559 

258 
290 
382 
547 
119 
125 
35 
54 

343 


47 

28 

19 

37 

76 

115 

129 

107 

53 

10 

-7 

384 

-191 

1146 

-269 

-590 

189 

539 

468 

285 

189 

100 

-61 

80 

91 

131 

231 

255 

250 

235 
278 
242 
603 
441 

272 
221 
322 
275 

236 
226 

273 
191 
183 
223 
288 


ERROR 

14 

36 

65 

16 

11 

1 

-10 

1 

109 

395 

96 

-299 

530 

-938 

349 

706 

-48 

-406 

-343 

-127 

-25 

-18 

264 

59 

42 

67 

-80 

015 

-75 

124 

-> 

1011 

-330 

-291 

-15 

338 

-64 

15 

146 

321 

-154 

-66 

-98 

-169 

t»5 
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TABLE  VI  (CONTINUED) 


WEEK 


46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 


LOSSES 

ESTIMATED  LOSSES 

ERROR 

118  ' 

435 

-317 

96 

453 

-357 

80 

411 

-331 

262 

339 

-77 

131 

340 

-209 

85 

330 

-245 

177 

342 

-165 

106 

389 

-283 

435 

387 

48 

145 

422 

-277 

91 

321 

-230 

126 

197 

-71 

442 

118 

324 

86 

178 

-92 

121 

180 

-59 

153 

220 

-67 

327 

261 

66 

191 

306 

-115 

220 

275 

-55 

336 

234 

102 

280 

237 

43 

237 

256 

-19 

202 

232 

-91 

300 

332 

-32 

273 

374 

101 

389 

366 

23 

646 

339 

307 

732 

331 

401 

359 

329 

30 

400 

268 

132 

253 

258 

-5 

247 

253 

-6 

249 

265 

-16 

248 

272 

-24 

536 

260 

276 

249 

293 

-44 

256 

265 

-9 

254 

269 

-15 

315 

312 

3 

208 

393 

-185 

L69 

448 

-179 

203 

484 

-281 

290 

464 

-174 

608 

419 

189 

357 

404 

-47 

260 

346 

-36 

189 

301 

-112 

166 

283 

-117 

262 

287 

-25 

216 

309 

-93 

224 

301 

-77 
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TABLE  VI  (CONTINUED) 


WEEK 


losses  ESTIMATED  LOSSES 


ERROR 


97 

304 

273 

98 

493 

257 

99 

234 

290 

100 

262 

308 

101 

466 

363 

102 

367 

459 

103 

357 

510 

104 

414 

513 

105 

393 

476 

106 

243 

406 

107 

438 

318 

108 

208 

295 

109 

404 

274 

110 

948 

310 

111 

735 

409 

112 

1056 

444 

113 

787 

485 

114 

1333 

463 

115 

453 

503 

116 

563 

451 

117 

475 

448 

118 

676 

459 

119 

582 

492 

31 

236 

-56 

-46 

103 

-92 

-153 

-99 

-83 

-163 

120 

-87 

130 

638 

326 

612 

302 

876 

-50 

112 

27 

217 

90 


Kean  Absolute  Error  =  166 
Mean  Error  =  11.33 

Error  Variance  =  6.28  X  104 


47 


TABLE  VII 


KALMAN  FILTER  RESULTS 
GROWING  SINE  WITH  HARMONICS 
CYCLE  LENGTH  32  WEEKS 


WEEK 

LOSSES 

ESTIMATED  LOSSES 

ERRCR 

1 

61 

60 

1 

2 

64 

63 

1 

3 

84 

63 

21 

4 

53 

61 

-8 

5 

87 

62 

25 

6 

116 

69 

47 

7 

119 

82 

37 

8 

108 

102 

6 

9 

167 

127 

40 

10 

405 

156 

249 

11 

*  89 

184 

-95 

12 

85 

-75 

160 

13 

339 

179 

160 

14 

208 

1155 

-947 

15 

80 

617 

-537 

16 

116 

969 

-853 

17 

141 

-362 

503 

18 

133 

162 

-29 

19 

125 

246 

-121 

20 

158 

183 

-25 

21 

164 

173 

-9 

22 

82 

117 

-35 

23 

203 

-102 

305 

24 

139 

215 

-76 

25 

133 

153 

-20 

26 

198 

165 

33 

27 

151 

345 

-194 

28 

140 

266 

-126 

29 

175 

227 

-52 

30 

359 

291 

68 

31 

280 

608 

-328 

32 

1253 

505 

748 

33 

273 

1782 

-1509 

34 

150 

1007 

-857 

35 

257 

258 

-1 

36 

559 

-12 

571 

37 

258 

198 

60 

38 

290 

32 

258 

39 

382 

21 

361 

40 

547 

144 

403 

41 

119 

384 

-265 

42 

525 

156 

-31 

43 

83 

50 

35 

44 

54 

-21 

75 

45 

343 

-67 

410 

48 


WEEK 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 


TABLE  VII  (CONTINUED) 


losses  ESTIMATED  LOSSES  error 


118 

96 

80 

262 

131 

85 
177 
106 
435 
145 

91 

126 

442 

86 
121 
153 
327 
191 
220 
336 
280 
237 
202 
300 
273 
389 
646 
732 
359 
400 

253 

247 
249 

248 
536 

249 
256 

254 
315 
208 
269 
203 
290 
608 
357 
260 
189 
166 
262 
216 


102 

75 

47 

24 

89 

85 

73 

105 

124 

240 

278 

308 

351 

466 

480 

492 

495 

518 

482 

436 

412 

367 

315 

263 

244 
227 

245 
324 
411 
396 
384 
332 
279 
229 
184 
206 
168 
142 
127 
134 
133 
156 
177 
222 
325 
371 
391 
391 
381 
379 


16 

21 

33 

238 

42 

0 

104 

1 

311 

-95 

-187 

-182 

91 

-380 

-359 

-339 

-168 

-327 

-262 

-100 

-132 

-130 

-113 

37 

29 

162 

401 

408 

-52 

4 

-131 

-85 

-30 

19 

352 

43 

88 

112 

188 

74 

136 

47 

113 

386 

32 

-111 

-202 

-225 

-119 

-163 
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TABLE  VII  (CONTINUED) 


WEEK  LOSSES  ESTIMATED  LOSSES  ERROR 


96 

224 

97 

304 

98 

493 

99 

234 

100 

262 

101 

466 

102 

367 

103 

357 

104 

414 

105 

393 

106 

243 

107 

438 

108 

208 

109 

404 

110 

948 

111 

735 

112 

1056 

113 

787 

114 

1339 

115 

453 

116 

563 

117 

475 

118 

676 

119 

582 

361 

-137 

340 

-36 

332 

161 

358 

-124 

334 

-72 

322 

144 

347 

20 

356 

1 

365 

49 

382 

11 

393 

-150 

374 

64 

384 

-176 

351 

53 

348 

600 

419 

316 

444 

612 

500 

287 

504 

835 

565 

-112 

504 

59 

460 

15 

414 

262 

404 

178 

Mean  Absolute  Error 
Mean  Error 
Error  Variance 


189.5 

7.64 

9.05  X  104 
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TABLE  VIII 


KALMAN  FILTER  RESULTS 
GROWING  SINE  WITH  HARMONICS 
CYCLE  LENGTH  52  WEEKS 


WEEK  LOSSES  ESTIMATED  LOSSES  ERROR 


1 

61 

63 

2 

64 

74 

3 

84 

82 

4 

53 

94 

5 

87 

-45 

6 

116 

-79 

7 

113 

14 

8 

108 

200 

9 

167 

272 

10 

405 

357 

11 

89 

738 

12 

85 

550 

13 

339 

90 

14 

208 

112 

15 

80 

-78 

16 

116 

-322 

17 

141 

-375 

18 

133 

-391 

19 

125 

-569 

20 

158 

-1997 

21 

164 

8194 

22 

82 

7445 

23 

203 

5585 

24 

139 

3774 

25 

133 

2474 

26 

198 

1713 

27 

151 

1335 

28 

140 

1058 

29 

175 

850 

30 

359 

706 

31 

280 

710 

32 

1253 

597 

33 

273 

1214 

34 

150 

799 

35 

257 

398 

36 

559 

192 

37 

258 

242 

38 

290 

65 

39 

382 

-47 

40 

547 

-90 

41 

119 

-75 

42 

125 

-295 

43 

85 

-537 

44 

54 

-989 

45 

343 

-3151 

-2 

-10 

2 

-41 

132 

195 

105 

-92 

-105 

48 

-649 

-465 

249 

96 

158 

438 

516 

524 

694 

2155 

-8030 

-7363 

-5382 

-3635 

-2341 

-1515 

-1184 

-918 

-675 

-347 

-430 

656 

-941 

-649 

-141 

367 

16 

225 

429 

637 

194 

420 

622 

1043 

3494 
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TABLE  VIII  (CONTINUED) 


WEEK 


LOSSES  ESTIMATED  LOSSES 


ERROR 


91 

92 

93 

94 

95 


118 

4022 

-3904 

96 

1478 

-1382 

80 

960 

-880 

262 

755 

-493 

131 

803 

-672 

85 

751 

-666 

177 

698 

-521 

106 

727 

-621 

435 

708 

-273 

145 

877 

-732 

91 

836 

-745 

126 

777 

-651 

442 

747 

-305 

86 

865 

-779 

121 

788 

-667 

153 

740 

-587 

327 

710 

-383 

191 

747 

-556 

220 

720 

-500 

336 

703 

-367 

230 

720 

-440 

237 

710 

-473 

202 

684 

-482 

300 

649 

-349 

273 

641 

-368 

389 

623 

-234 

646 

631 

15 

732 

689 

43 

359 

750 

-391 

400 

724 

-324 

253 

707 

-454 

247 

663 

-416 

249 

624 

-375 

248 

590 

-342 

536 

561 

-25 

249 

583 

-334 

256 

559 

-303 

254 

540 

-286 

315 

525 

-210 

208 

522 

-314 

269 

504 

-235 

203 

498 

-295 

290 

481 

-191 

608 

479 

129 

357 

527 

-170 

260 

527 

-267 

189 

509 

-320 

166 

478 

-312 

262 

443 

-181 

216 

423 

-207 
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TABLE  VIII  (CONTINUED) 


WEEK 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 
119 


LOSSES 


ESTIMATED  LOSS E$ 


224 

304 

493 

234 

262 

466 

367 

357 

414 

393 

243 

438 

208 

404 

948 

735 

1056 

787 

1339 

453 

563 

475 

676 

582 


395 

368 

351 

359 

328 

302 

303 
291 
280 
278 
276 
257 

270 
255 

271 
358 
414 
506 
560 
675 
677 
690 
689 
710 


ERROR 


-171 

-64 

142 

-125 

-66 

164 

64 

66 

134 

115 

-33 

181 

-62 

149 

677 

377 

642 

281 

779 

-222 

-114 

-215 

-13 

-128 


Nean  Absolute  Error  =  551 
Mean  Error  =  -359.16 

Error  Variance  =  17.47  X  105 
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TABLE  IX 


MEEK 

1 


KALMAN  FILTER  RESULTS 
LINEAR  MODEL 


LOSSES  ESTIMATED  LOSSES  ERROR 


61  -21  82 

64  -6  70 

84  10  74 

53  26  27 

87  41  46 

116  57  59 

119  73  46 

108  88  20 

167  104  63 

405  120  285 

89  135  -46 

85  122  -37 

339  127  212 

208  185  23 

80  209  -129 

116  203  -87 

141  203  -62 

133  206  -73 

125  206  -81 

158  204  -46 

164  208  -44 

82  210  -128 

203  200  3 

139  209  -70 

133  207  -74 

198  203  -5 

151  210  -59 

140  209  -69 

175  206  -31 

359  208  151 

280  234  46 

1253  247  1006 

273  286  -13 

150  294  -144 

257  298  -41 

559  305  254 

258  324  -66 

290  330  -40 

382  337  45 

547  348  199 

119  366  -247 

125  364  -239 

85  361  -276 

54  357  -303 

343  351  -8 
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TABLE  IX  (CONTINUED) 


TABLE  IX  (CONTINUED) 


WEEK 

LOSSES 

ESTIMATED  LOSSES 

ERROR 

96 

224 

374 

-150 

97 

304 

372 

-68 

98 

493 

372 

121 

99 

234 

380 

-146 

100 

262 

377 

-115 

101 

466 

376 

90 

102 

367 

382 

-15 

103 

357 

385 

-28 

104 

414 

387 

27 

105 

393 

391 

2 

106 

243 

394 

-151 

107 

438 

391 

47 

108 

208 

396 

-188 

109 

404 

392 

12 

110 

948 

395 

553 

111 

735 

416 

319 

112 

1056 

429 

627 

113 

787 

450 

337 

114 

1339 

462 

877 

115 

453 

485 

-32 

116 

563 

488 

75 

117 

475 

494 

-19 

118 

676 

497 

179 

119 

582 

505 

77 

Mean  Absolute  Error  =  141 


Mean  Error  =  -15.33 

Error  Variance  =  4.18  X  104 
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table  X 


WEEK 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 


KALMAN  FILTER  RESULTS 
CONSTANT  MODEL 


LOSSES  estimated  losses 


ERROR 


61 

64 

84 
53 
87 

116 

119 

108 

167 

405 

89 

85 
339 
208 

80 

116 

141 

133 

125 

158 

164 

82 

203 

139 
133 
198 
151 

140 
175 
359 
280 

1253 

273 

150 

257 

559 

25S 

290 

382 

547 

119 

125 

85 

54 

343 


29 

77 

80 

103 

103 

163 

164 
117 

94 

136 
208 
208 
149 
142 

137 
131 

157 

158 
153 
201 
198 
165 
198 
196 
171 

175 

176 
260 
267 
273 
264 

258 
271 

259 
287 
305 
317 
247 
164 
112 

7 


61 

64 

84 

52 

58 

39 

39 

5 

64 

242 

-75 

-32 

245 

72 

-128 

-92 

-8 

-9 

-12 

27 

7 

-76 

50 

-62 

-65 

33 

-47 

-56 

4 

184 

104 

993 

6 

-123 

-7 

301 

-13 

31 

?5 

242 

-198 

-122 

-79 

-58 

266 
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TABLE  X  (CONTINUED) 


WEEK 

LOSSES 

ESTIMATED  LOSSES 

ERROR 

46 

118 

123 

-5 

47 

96 

118 

-22 

48 

80 

no 

-30 

49 

62 

96 

166 

50 

131 

141 

-10 

51  i 

85 

133 

-48 

52  ’ 

177 

120 

57 

53 

1C6 

176 

-70 

54 

435 

174 

261 

55 

145 

228 

-83 

56 

91 

166 

-75 

57 

126 

125 

1 

58 

442 

126 

316 

59 

86 

132 

-46 

60 

121 

90 

31 

61 

153 

120 

33 

62 

327 

140 

187 

63 

191 

141 

50 

64 

220 

191 

29 

65 

336 

194 

142 

66 

280 

198 

82 

67 

237 

270 

-33 

68 

202 

264 

-62 

69 

300 

242 

58 

70 

273 

299 

-26 

71 

389 

298 

91 

72 

646 

360 

286 

73 

732 

371 

361 

74 

359 

505 

-146 

75 

400 

401 

-1 

76 

253 

400 

-147 

77 

247 

391 

-144 

78 

249 

317 

-68 

79 

248 

266 

-18 

80 

536 

252 

284 

81 

249 

268 

-19 

82 

256 

250 

6 

83 

254 

255 

-1 

84 

315 

254 

61 

85 

208 

305 

-97 

86 

269 

304 

-35 

87 

203 

277 

-74 

88 

290 

254 

36 

89 

608 

290 

318 

90 

357 

290 

67 

91 

260 

356 

-96 

92 

189 

353 

-164 

93 

166 

302 

-136 

94 

262 

223 

39 

95 

216 

262 

-46 
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TABLE  X  (CONTINUED 


WEEK 


LOSSES  estimated  losses 


96 

97 
96 
99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 
119 


224 

304 

493 

234 

262 

466 

367 

357 

414 

393 

243 

438 

208 

404 

948 

735 

1056 

787 

1339 

453 

563 

475 

676 

582 


262 

241 

296 

302 

254 

261 

267 

355 

356 
400 
398 
391 
438 
437 
406 
415 
660 
813 
788 
798 
581 
563 
548 
637 


Mean  Absolute  Error  = 

109 

Mean  Error  = 

41.33 

T  ***or  V?— 

Q  ?3  X  103 

ERROR 


-38 

63 

197 

-68 

8 

205 

100 

2 

58 

-7 

-155 

47 

-230 

-33 

542 

320 

396 

-26 

551 

-345 

-18 

-88 

128 

-55 
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TABLE  XI 


LINEAR  REGRESSION  RESULTS 


WEEK 

LOSSES 

ESTIMATED  LOSSES 

ERROR 

13 

339 

196 

143 

14 

208 

252 

-44 

15 

80 

256 

-176 

16 

116 

224 

-108 

17 

141 

208 

-67 

18 

133 

201 

-68 

19 

125 

193 

-68 

20 

158 

184 

-26 

21 

164 

184 

-20 

22 

82 

185 

-103 

23 

203 

170 

33 

24 

139 

179 

-40 

25 

133 

176 

-43 

26 

198 

172 

26 

27 

151 

178 

-27 

23 

140 

177 

-37 

29 

175 

174 

1 

30 

359 

177 

182 

31 

280 

203 

77 

32 

1253 

217 

1036 

33 

273 

350 

-77 

34 

150 

351 

-201 

35 

257 

337 

-80 

36 

559 

336 

223 

37 

258 

369 

-111 

38 

290 

366 

-76 

39 

382 

366 

16 

40 

547 

376 

171 

41 

119 

402 

-283 

42 

125 

383 

-258 

43 

85 

366 

-281 

44 

54 

347 

-293 

45 

343 

327 

16 

46 

118 

334 

-216 

47 

96 

320 

-224 

48 

80 

306 

-226 

49 

262 

291 

-29 

50 

131 

292 

-161 

51 

85 

283 

-198 

52 

177 

271 

-94 

53 

106 

266 

-160 

54 

435 

257 

178 

55 

145 

272 

-127 

56 

91 

265 

-174 

57 

126 

255 

-129 
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TABLE  XI  (CONTINUED) 


WEEK 


LOSSES  ESTIMATED  LOSSES 


ERROR 


58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 
103 


442 

248 

194 

86 

263 

-177 

121 

253 

-132 

153 

24b 

-93 

327 

241 

86 

191 

248 

-57 

220 

246 

-26 

336 

246 

90 

280 

253 

27 

237 

256 

-19 

202 

256 

-54 

*’00 

255 

45 

J73 

259 

14 

389 

261 

128 

646 

270 

376 

732 

293 

439 

359 

319 

40 

400 

323 

77 

253 

330 

-77 

247 

329 

-82 

249 

327 

-78 

248 

326 

-78 

536 

324 

212 

249 

337 

-88 

256 

336 

-80 

254 

334 

-80 

315 

333 

-18 

208 

334 

-126 

269 

331 

-62 

203 

330 

-127 

290 

327  ' 

-37 

608 

327 

281 

357 

342 

15 

260 

345 

-85 

189 

343 

-154 

166 

339 

-173 

262 

334 

-72 

216 

333 

-117 

224 

330 

-106 

304 

327 

-23 

493 

328 

165 

234 

337 

-103 

262 

334 

-72 

466 

333 

133 

367 

340 

27 

357 

343 

14 
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TABLE  XI  (CONTINUED) 


WEEK 

LOSSES 

ESTIMATED  LOSSES 

ERROR 

104 

414 

346 

68 

105 

3S3 

350 

43 

106 

243 

354 

-111 

107 

438 

352 

86 

108 

208 

357 

-149 

109 

404 

353 

51 

110 

948 

357 

591 

111 

735 

380 

355 

112 

1056 

395 

661 

113 

787 

421 

366 

114 

1339 

437 

902 

115 

453 

472 

-19 

118 

563 

474 

89 

117 

4/b 

480 

-5 

118 

676 

484 

192 

119 

582 

493 

89 

Mean  Absolute  Error  =  139 


Mean  Error 
Frror  Variance 


*  9.9 
4.53  X  104 
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2 
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120 


)  FORECAST  MEEKS 


V.  CONCLUSIONS 


In  the  case  of  the  linear  model.  It  was  shown  In  Chapter  III  that 
E[C(t)]  -  (ai+a2t)fT(l)[I-LT+h  fT(l)]_1  h 

-  a2  fT(l)[I-LT+h  f^nrV-h  fT(l)][I-L7+h  fT(l )3_1  h 

for  large  t  and  when  the  actual  data  was  such  that 
E[C(t)]  =  a1  +  a2t. 

Then,  E[C(t)]  =  a^  +  a2t  +  a2* 

This  showed  that  when  the  assumed  model  was  correct,  the  forecast  was 
biased  by  the  slope  of  the  data  line.  In  the  case  of  the  problem  of 
forecasting  the  Marine  losses,  the  least  squares  line  through  the  one 
hundred  and  nineteen  points  had  slope  3.41.  This  suggested  that  Exponen- 
1 1  a  1  should  have  produced  a  forecast  which  Is  3.41  more  than  the  actual 
data,  or  an  error  of  -3.41.  The  actual  mean  error  was  -2.21.  Also,  the 
Exponential  Smoothing  forecast  would  be  above  the  linear  regression  fore¬ 
cast  when  the  losses  were  Increasing  and  below  when  they  were  decreasing. 
The  results  bear  this  out,  although  not  around  the  time  periods  when  the 
actual  losses  changed  directions.  For  example,  the  losses  Increased 
during  the  time  periods  thirty-seven  through  forty.  The  Exponential 
Smoothing  forecast  with  the  linear  model  exceeded  the  linear  regression 
forecasts  during  these  periods.  The  losses  then  decreased  during  the 
periods  forty-one  through  forty-four.  The  Exponential  Smoothing  forecast 
was  above  the  linear  regression  forecast  for  the  forty-first  week  but  It 
then  was  below  the  linear  regression  forecasts  through  forty-four. 
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Standard  ,nn 
Deviation  JUU 


ESI  6 
4 

+  4 
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KF 

a 
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The  results  obtained  with  the  Kalman  Filter  and  the  growing  sine  model 
were  generally  comparable  with  the  results  with  Exponential  Smoothing. 

This  Is  evidenced  in  Figure  7.  The  exception  was  the  sine  model  with  a 
cycle  length  of  fifty- two  weeks.  The  Kalman  Filter  forecast  for  this 
model  fluctuated  badly  for  the  Initial  forecasts.  However,  the  last 
seventy  weeks  of  the  data  produced  a  mean  absolute  error  of  -321.71  and 
a  mean  error  of  -208.71.  The  last  forty  weeks  produced  a  mean  absolute 
error  of  219  and  a  mean  error  -24.08.  When  these  results  are  compared 
with  the  overall  mean  absolute  error  of  651  and  the  mean  error  of  -359.16, 

It  appears  that  the  forecasts  were  Improving  as  the  effects  of  the  Initial 
forecasts  were  dampening  out.  Thus,  the  model  whould  not  be  totally  dis¬ 
missed,  especially  since  the  Exponential  Smoothing  results  with  the  cycle 
length  of  fifty- two  weeks  did  not  vary  greatly  from  the  results  when  a 
cycle  length  of  sixteen  and  of  thirty-two  weeks  was  used.  The  Kalman 
Filter  technique  assumes  that  the  forecast  error  variance  Is  known.  For 
the  problem  Investigated,  that  variance  was  estimated,  which  was  a  variation 
on  the  Kalman  Filter.  If  the  true  variance  had  been  known,  the  results  of 
all  the  models  with  which  the  Kalman  Filter  was  used  might  have  been  Im¬ 
proved. 

With  the  models  used  In  forecasting  Marine  losses,  the  simpler  the 
model,  the  better  the  results.  Captain  O'Brien  found  that  the  linear  model 
produced  better  results  than  the  more  complicated  growing  sine  model  for 

the  prediction  of  monthly  Marine  losses  [O'Brien  1968].  The  same  was  true 

t> 

for  the  weekly  forecasts.  In  addition,  the  even  simpler  constant  model 
produced  even  better  results.  However,  If  the  actual  cycle  period  had  been 
used  rather  than  an  estimate,  the  results  of  the  Exponential  Smoothing  with 
the  growing  sine  model  might  have  been  greatly  Improved.  The  period  length 
can  be  determined  as  described  In  R.G.  Brown's  book[Brown  1963,  p.  66-75]. 


In  general,  the  estimated  mean  of  the  forecast  errors  was  not  zero. 
For  example,  the  Kalman  Filter  with  the  constant  model,  on  the  average, 
forecasted  forty-one  fewer  losses  than  actually  occurred.  With  this 
Information,  a  better  policy  would  have  been 

C'(t)  =  C(t)  +  u, 

where  u  Is  the  mean  forecast  error,  C'(t)  Is  the  estimated  losses  and 

A 

C(t)  Is  the  forecasted  losses  from  the  method  used.  The  variance  of  the 
forecast  error  will  not  be  changed  by  this  modification  since  C(t)  would 

a 

still  be  computed  using  C(t-1)  -  C(t-l).  However,  the  mean  error,  where 
the  error  Is 

e  * (t)  *  C(t)-C'(t), 
will  be  zero. 

Finally,  the  results  Indicated,  as  evidenced  In  Figure  7,  that  even 
when  an  estimate  of  the  forecast  error  covariance  matrix  Cy  Is  used,  the 
Kalman  Filter  results  do  not  differ  greatly  from  those  of  Exponential 
Smoothing.  Since,  except  for  the  constant  model,  the  Kalman  Filter  Is 
more  complicated  than  either  the  Exponential  Smoothing  or  the  Linear 
Regression  methods,  either  of  these  latter  methods  should  be  used  when 
Cy  Is  unknown  or  when  no  computer  facilities  are  available.  When  Cy  Is 
known  and  a  computer  Is  available,  the  Kalman  Filter  technique  might 
provide  better  results  than  either  of  the  other  two  methods. 

Under  the  criterion  chosen,  the  minimum  mean  absolute  error,  the 
Kalman  Filter  method  using  the  constant  model  was  the  best  forecasting 
method.  This  method  could  be  easily  used  for  forecasting  Marine  losses 
In  Vietnam.  The  calculations  for  this  constant  model  are  relatively 
simple  and  can  be  accomplished  with  a  calculator  without  the  aid  of  a 
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computer.  The  calculation  of  the  estimate  of  the  forecast  error  variance 
Is  the  most  complicated  portion  of  the  constant  model.  The  mean  absolute 
error  and  the  variance  of  the  forecast  errors  achieved  with  this  model 
compare  favorably  with  the  data  mean  of  294.6  and  variance  of  the  data  of 

4 

5.21  X  10  .  Using  this  model  would  be  more  beneficial  than  merely  using 
the  data  mean  as  the  forecast  of  the  losses  for  the  next  time  period. 
Because  of  the  difference  between  the  two  variances,  one  would  be  more 
certain  of  a  minimum  error  with  the  Kalman  Filter  constant  model. 


GENERAL  EXPONENTIAL  SMOOTHING  PROGRAM 
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This  program  computes  successive  weekly  forecasts  for  the  model 
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•Mil 


DATA. e/e. 86603. 0,94868, 0.97468/ 


ft  1 
1$: 

40(5,6)  (X(II, 1-1,119) 


MMM}, 1-1,3} 


LW-11 . 

DO  91  J-1,119 
91  TIME! Jl-J 

So  1  ftl;i 

MMM 
R 

1  Fl 

Rc-UISfOf  IXIi 
6  FORMAT ( 11F6.0 ) 

HjtfikhS 

?fL|  -MS1*""" 

r]l!  -o.o 

21  CONTINUE 
00  8  1-1,3 
ESUM-0.0  , 

CALCULATE  ESTIMATE  ANO  ERROR  OF  ESTIMATE 

JO  lg  IT-l.LM 

DC  40  LK-l.MMM 
5*S-»A£LKI*rONEILK ) 

40  CONTINUE 

E?IT j-Xt if l-XEST ( IT  I 

mmsw* 

RELERf IT}  «  EIITI/XESTIIT} 

UPCATI  COEFFICIENTS 

. 

R(M I  ♦  UCMI 


••II,*  WEIGHTING  CONSTANT.  8  «  •, 


74 


XREAN*  SIGMA* DEV) 


00  202  NN«1*LN 

202 

ABERR-SUMAi/LM 

losses'-10*- 

>»  gffilSF*- 


ERROR 


SF4.0) 


,S2*v§55ai*ii!vS,?;gCA*  error  »  •*Fi.2,//.iox, 

*  '  F?0  5|^ANCE  *  SE14.7,//,  10X*  *StInOaIo  OBVIATION 

UrItEi A* 201) 

20l*LCSSFj.<IJh?8,,,'PL0T  CF  TI"E  VS  L0SSES  RNO  ESTIMATED 

Wkmkvini^ 

SRITEU?20|)IME,X,LP,3, 

*°j  g^2^)*iiHi:tta.K eMms  **•««•••///» 

177  pR||giiF0We*2**1-0 


•FONE<2)41.0 


DO  io  I-l.L 
1°  A  -  AjXltl 

„  DC  20  I-l.L 
20  I  *  gWxflf-A) 

LiufcfeT*B’ 

END 


KALMAN  FILTER  PROGRAM 


This  program  computes  successive  weekly  forecasts  for  the  growing 
sinusoidal  model  with  a  cycle  length  of  32  weeks. 
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This  paper  Investigates  the  application  of  the  Kalman  Filter  and  the  General 
Exponential  Smoothing  techniques  of  forecasting.  Both  methods  are  derived  and 
the  similarities  and  differences  between  them  are  discussed.  The  two  techniques 
are  then  applied  to  the  practical  problem  of  predicting  weekly  losses  suffered 
by  the  U.  S.  Marine  Corps  units  in  the  I  Corps  Tactical  Zone  In  the  Republic 
of  Vietnam.  The  mean  absolute  error  of  the  prediction  Is  used  as  the  criterion 
for  shooslng  the  better  of  the  two  methods.  Results  are  given  for  both  techniques 
as  well  as  for  the  method  of  linear  regression.  In  general  the  Kalman  Filter 
provides  the  smallest  mean  absolute  error  for  the  three  mathematical  models; 
linear,  growing  sine  with  harmonics  and  frequency  of  sixteen,  thirty-two,  and 
fifty-two  weeks,  and  a  constant  model. 
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